Full Counting Statistics for Number of Electrons Dwelling in a Quantum Dot 
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Motivated by recent real-time electron counting experiments, we evaluate the full counting statis- 
tics (PCS) for the probability distribution of the electron number inside a quantum dot which is 
weakly coupled to source and drain leads. A non-Gaussian exponential distribution appears when 
there is no dot state close to the lead chemical potentials. We propose the measurement of the joint 
probability distribution of current and electron number, which reveals correlations between the two 
observables. We also show that for increasing strength of tunneling, the quantum fluctuations qual- 
itatively change the probability distribution of the electron number. In this paper, we derive the 
cumulant generating functions (CGFs) of the joint probability distribution for several cases. The 
Keldysh generating functional approach is adopted to obtain the CGFs for the resonant-level model 
and for the single-electron transistor in the intermediate conductance regime. The general form for 
the CGF of the joint probability distribution is provided within the Markov approximation in an 
extension of the master equation approach [D. A. Bagrets, and Yu. V. Nazarov, Phys. Rev. B 67, 
085316 (2003)]. 

PACS numbers: 73.23.Hk,72.70.-|-m 
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I. INTRODUCTION 

Measurements of the average current and its fluc- 
tuation (noise) have been powerful tools to study the 
quantum transport in mesoscopic systems^. In the last 
decade, the theory of full counting statistics (FCS)^ has 
been established. The FCS provides a probability distri- 
bution of the current P{I) from which one obtains not 
only first and second cumulants, i.e. the average current 
and the noise, but also any order cumulant. Precisely, 
/ = q/to is the time-averaged current during the mea- 
surement time to, where 



to/2 



dt I{t)/e, 



(1) 



-to/2 



is the number of transmitted electrons. 

The last few years, there have been remarkable ad- 
vances in the experimental study of the FCS^'^'^'^. The 
third cumulant of the probability distribution of current 
has been measured for a tunnel junction'^-^. Recently, 
the probability distribution itself was experimentally ob- 
tained for a quantum dot (QD) coupled to source (R) and 
drain (L) leads^'^. The experiments adopted the real- 
time electron-counting technique using quantum point 
contact (QPC) charge detectors. The measured quan- 
tity is the time evolution of electron number inside the 
QD, n{t). The experiment*^ demonstrated that n{t) fluc- 
tuates between and 1, as schematically shown in Fig.[T] 
The number of transmitted electrons through the QD, q, 
is obtained by counting the number of transitions from 
to 1 [kinks pointed by arrows in Fig.[T]. However, the 
technique works when electrons tunnel in one direction 
only, in which case there is a one-to-one correspondence 
between the changes of n{t) from to 1 and the tunnel- 
ing processes from one lead L into the QD. Therefore, 
the analysis of the experiment^ has been limited to large 



source-drain voltages T <C eV (we use fee = ^ = 1) where 
the thermal fluctuations are suppressed. 

The experiments were also performed in the weak tun- 
neling case, where the FCS theory within the Markov ap- 
proximatioiiiil:-^ works perfectly. For increasing strength 
of the tunnel coupling, the FCS is particularly interest- 
ing because quantum coherence effects show up. How- 
ever, the method used in the experimenijfi' does not work 
for this case either, since electrons tunnel back and forth 
coherently. Therefore it is important to extend the FCS 
theory and to relate it directly with quantities that are 
measurable in the QPC experiments. In this paper, we 
will provide the FCS for the time-averaged electron num- 
ber N = r/io, where 



to/2 

dt n(t), 

-to/2 



(2) 



can be obtained from the experimental data. The observ- 
able N is well defined for any condition including small 
source-drain voltages T ^ eV. On top of that, from 
the experimental setup in Ref.^, one can obtain more in- 
formation, namely the joint probability distribution of 
current and electron number P{I,N), since at a large 
source-drain voltage T<$^eV, the data from a single mea- 
surement [Fig. [1] provides particular values for both of q 
and r. 

In what follows, we will calculate the characteristic 
function Z or the cumulant generating function (CGF) 
W for the joint probability distributions, 

/oc 
drP(/,Y)e^«^ + "«, (3) 

'^''^'^'^hsqHr'y (4) 



W(A,0-lnZ(A,0=^ 
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probability distribution of current only by Bagrets and 
NazarovM. 

The structure of the paper is the following. In Sec. [TTl 
we discuss the probability distribution of electron num- 
ber for the spinless non-interacting QD, the RLM. Then 
in Sec. IIIIl the spin effect in an interacting QD is con- 
sidered within the Markov approximation. In Sec. II VI we 
present the joint probability distribution of current and 
electron number. In Sec. |V1 we will demonstrate how 
the quantum fluctuations of charge affect on the statis- 
tical properties of electron number inside the QD. We 
summarize in Sec. I VII 



FIG. 1: The time evolution of electron number for rare (upper 
panel) and frequent (lower panel) tunneling cases. 



The parameters A and ^ are called the 'counting fields'. 
The aim of this paper is to analyze basic features of the 
statistical properties of electron number inside the QD 
and the joint probability distribution for several relevant 
cases. 

In order to derive the CGFs for the joint probabil- 
ity distribution, we adopt two approaches, the Keldysh 
generating functionat^ i^°'^^ and the master equation ap- 
proachesii^. The former has been used widely in meso- 
scopic physicsi^ii^ii^ii^ii^iiiii^ii^i2£i2ii2^ independently of 
the FCS study. However, they are closely related to 
each other: The Keldysh generating functional is reduced 
to the CGF by fixing the 'quantum' components of the 
source fields as the counting fieldsiS,. In sections [II Al and 
IV Al we extend our previous Keldysh generating func- 
tional study for the current and charge noises^^ to the 
FCS. In the formal level, our extension could be consis- 
tent with the FCS theory of charge developed by Pil- 
gram and Biittiker— . Their scheme is based on Levi- 
tov's 5'-matrix approach- and suitable to deal with an 
open and noninteracting chaotic cavity, since it provides 
a systematic method for the disorder averaging^!. At 
the same time, the scheme requires more involved calcu- 
lations to account for the effects of the back scattering 
from the tunnel barriers and the Coulomb interaction. 
Therefore ReL^^ analyzes the situations where both of 
them are weak. The present paper will address comple- 
mentary cases, closed and interacting QDs. In Sec. Ill Al 
we consider a simple model for a non-interacting QD, 
the resonant-level model (RLM). Through the calcula- 
tions we would like to demonstrate that the Keldysh 
generating functional is a useful tool for the FCS study. 
In Sec. IV Al we consider an interacting QD, the single- 
electrons transistor (SET) in an intermediate conduc- 
tance regime. There the approximate generating func- 
tional derived in Refsi ^°'^^ is extended to the CGF. In 
Sec. nil Al we also derive the general CGF for the joint 
probability distribution of current and electron num- 
ber within the Markov approximation: We extend the 
master equation approach originally developed for the 



II. FCS FOR A NON-INTERACTING QD 

In this section, we derive the CGF for the RLM. Our 
calculational tool is the Keldysh generating functional in 
the path-integral representation (We refer to Refsj ^°'^^ 
for details). We will check that our result is consis- 
tent with previous theories up to the second cumulant 
(Sec. Ill A| . In Sec. Ill B[ we discuss the equilibrium prob- 
ability distribution of electron number in the weak tun- 
neling regime. 



A. FCS formalism based on the Keldysh 
generating functional 

The Hamiltonian of a non-interacting closed QD for 
spinless case (RLM) is given by, 

H = ^ ekalf^hrk + eo d^ d + ^{Vrd^cirk + H.c.) , (5) 



where a^fc annihilates an electron with wave vector k in 
the lead r while, d annihilates an electron in the QD 
with the dot- level eo- The electrons in the lead r obey a 
Fermi distribution /(w - ^lr) = l/(e('^-^'-)/^ + 1). The 
chemical potential /z^ is fixed by the source-drain voltage 

In order to derive the CGF for the joint probability dis- 
tribution of q and t, equivalently / and iV, we first intro- 
duce the closed time-path C (Fig. [2]). The choice of this 
time-path is particularly convenient since with the help of 
the technical know-how to perform the path integra l "'^°i^^ 
and to solve the differential equation for the Green func- 
tion (GF) defined on C^, such as Eq. p^ . the calcula- 
tions of the characteristic function become quite similar 
to those of the partition function in the imaginary-time 
formalism. Next, we introduce two source fields, Lpr{t) 
and h{t), the phase of the tunneling matrix element and 
the fiuctuation of the dot-leveP^: 



Vr 



y *V.(t) 



eo 



eo - h{t). 



(6) 



Note that here the time t is defined on the closed time- 
path C. After project the time on C onto the real axis, 



two components, fr+{h^) and Lpr-{h-) residing Cj^ and 
C_, appear. The averages of the two components, 

ifr+{t)+Vr-{t) _ , h+(t)+hjt)_ 

posses the physical meaning. The differences are ficti- 
tious degrees of freedom called the 'quantum' compo- 
nentsiS,. In the generating functional approach, these 
components are set to in the end of calculations^ i^°i^^ . 
In the FCS scheme, they are fixed as constants, the count- 
ing fields. 



^Pr+{t)- ipr-{t) =\r, h+{t) - h^{t) = S,, 



(8) 



during the measurement (— to/2 <t < to/2) and other- 
wise. 

We introduce the path-integral representation of the 
Keldysh generating functional from the Hamiltonian ([5]) 
by the time-slicing metho d^°i^^ . Under the condition ([8]) , 
it is reduced to the CGF: 



yV[\i) = \n V[a;^,,ark,d*,d] exp (iS") , 



(9) 



^ ^ dt{y^ark{t)*{idt~ £k)ark{t) 

-fc y rk 

+ d{t)*[idt-ea + h{t)]d{t) 

- ^(V^.e^'^-(*)d(t)*a,.fc(t)+c.c.)l, 

rk ) 



(10) 



where Grassmann fields satisfy the antiperiodic boundary 
condition, for example d[— oo— 1(3) = — d{^oo). Although 
the counting fields \l and A^ appear at this stage, the 
final result contains only A = Al — Xn from the charge 
conservation^*'. We can check the first derivative in ^ 
gives r with an offset; 



dW 



dm 



A=4=0 



'" dt([dt(t),rf(t)])=r-|, (11) 

-to/2 ^ 



where d{t) is in the Heisenberg picture and the average 
is performed over the Hamiltonian ([5]) without the tun- 
neling term (For technical details, see Appendix. A of 
Refi^). From now on, we will shift N by the offset —1/2. 



c+ 

-^ — 



^'Ct 



-iy3 



FIG. 2; The closed time-path consisting of forward C+, back- 
ward C- and imaginary time Ct branches. 

After the path integral over ark along GlSiii, which is 
formally to complete the square, we obtain the effective 



action for the QD as (In the following, we will omit trivial 
constants) , 

Sd= f dtdt' d{t)* [g-^{t,t') ~T.^^^{t,t')]d{t'). (12) 
Jc 

The GF satisfies the following equation and the anti- 
periodic boundary condition: 



{idt-eo)9{t,t')^5{t,t'), 
.g(— oo + ij3, t') = —g{~(X),t'). 



(13) 
(14) 



The (5-function is defined on C and satisfies 
J^dt6{t,t') = 1. Precisely g is a discrete matrix, 
since time is discreteiSiiii. In the continuous notation, 
we can regard Eq. p3)) as the differential equation 
defined on C, which is given as a function t = z{s) with 
s, a real and monotonically increasing parameter. Then 
Eq. (fT51) can be solved^SiS^ by considering the anzata^^, 

g{t, t') = g-+{t - t') e{t, t') + g+-{t - t') e{t', t), (15) 

where the step function is defined on C: dt9{t,t') = 
—df 0{t,t') = S{t,t'). We observe that g^^ becomes. 



^-+ 



The anti-periodic boundary condition p4p . g '"(— oo — 
iP,t') = —g^ (— oo,i'), determines the parameter a as 
a = /(eo). After projecting the time on C onto the real 
axis, Eq. fTS]) is mapped onto the 2x2 Keldysh space: 

'g++{t-t') g+-{t~t') 
_g-+{t-t') g-{t~t') 

+- (t-t')0{t-t')+g-+{t-t')e{t' -t), 
+-{t-t')0{t'-t)+g-+{t-t')0{t-t'). 

A compact representation follows from the Keldysh ro- 
tation by an orthogonal matrix Q = (t° — «t^)/-\/2 as 
g = QgQ~^ [ "J"* (* = li2, 3) are the Pauli matrices 
and T*' is the unit matrix ]. Then in Fourier space, 
the matrix GF is written with the retarded, g^{uj) = 
l/{iLj + iO — Co), the advanced g^ = g^* and the Keldysh 
g^{uj) = 2iIm5^(Lj) tanh(a;/2r) components as. 





g(t, 


t') = 


g++(t- 
g—{t- 


-t') = g 
-t') = g 



g^ioj) 
g^icu) g^(c.) 



gi^) 

The self-energy in Eq. p2|) contains the source fields: 
E^^«(t,t') =J2 S^(t,t')e'^*''-^*^~^'-^*'»-/i(t)<^(t,t')- 

r=L,R 

The 'bare' part Er{t,t') =X]fc I KfcPsr-fc (*,*') is written 
with the lead-electron GF g^k, which is obtained from 
g by replacing eq with Srk- The retarded and Keldysh 
components of S^ , 






-iTr tanh ■ 



2T' 



(16) 
(17) 



are then written with the couphng strength F^ = 
27r^j, lyrpf^Cw — Sk), which is assumed to be energy in- 
dependent. 

Further performing the path integral over d, we obtain, 



where F^ 



W = Tr In 



(zG^-«"'), G^-«"'=5-i-E^-«, (18) 



where the trace is understood as the time integrations 
along C. It can be performed utilizing the asymptotic 
form, 



Jc 
2 



dtTi 

to/2 



G^Ht,t) 



!£Tr 
2 



G^'«(0) 



for to^oo. After the Fourier transform and the integra- 
tion in terms of ^, we obtain, 



G^-Hu 



W = —f du;TrliiG^-^{Lu)-\ 
2ttJ 

y^ = {uj - eo)Ti+ -To 

_ ^ e''^'-^i/ViE,(cc;)Tie-'^'-'^i/2. 



r=L,R 



The frequency in S^l^) is shifted as w — > w — /x^ in 
order to absorb the part Urt in the phase, Eq. ([7]). After 
straightforward calculations and subtracting a constant 
in order to fulfill the normalization condition W(0,0) = 
0, we come to the final form: 



^ r=L,R ' ■' 

T{uj) = A\G"iLu)flm^f{uj)lmJ:^{uj), 
G^iu;) - 



1 



-eO-Er-=L,fl^?('^) 



(19) 
(20) 

(21) 



where the transmission probability T{uj) ~ Ti,Tii/[(uj — 
eo)^ +rV4] is the Lorentzian (F = Fl +Ffl,). The fourth 
line of Eq. ITOl) is our new result. Without the fourth 
line, Eq. (fTOl) is reduced to the Levitov-Lesovik formula^. 
Further, the integration in frequency can be performed 
for T :^ r following the procedure in Sec. IV of Ref4, 



w'^'Hx,0=tor: 



y/D~l 



D = 1 



,F+F 



L^ fl /AX 



F2 



(e'^-l) + 4 



Fs=F^ 

^ _»> L 

F2 



(22) 



(c— -1) 






Z-^r 



=L,R 



r^ and 



F±=F,/(±eoT/ir), 



(23) 



is the tunneling rate of an electron into/out of the QD 
through the junction r within Fermi's golden rule. Note 
that the condition T^F is that for the Markov approx- 
imation and that Eq. ((22)) is the first order expansion of 
Eq. (fT9|) in F assuming ^ ex F. Actually, the first and 
second cumulants, {{St)) and {{St'^)), and the covariance 
{{SqSr)), 



{{Sr))/to = 
{{Sr'))/h = 



{{SqSr)) 



F+-F- 
2Fs ' 
F+F" 

, {{Sq)){{Sr)) 
Fs to ' 



(24) 
(25) 
(26) 



are reproduced by the master equation approach for the 
noise^i. Moreover, later in Sec. IIII Ai we will rederive 
Eq. (|22p in an extension of the master equation ap- 
proach. The second check is to compare the expres- 
sions for the second cumulants and the covariance ob- 
tained from Eq. (|19p with those obtained by the stan- 
dard Keldysh diagrammatic approach for the Anderson 
model2^. We check our results, 



{{Sr'))/to 
{{STSq))/to = F 



p2 I duj /cff(c^)[l-/cff(^)] ^27) 



27r[(c^-eo)2 + F2/4]2' 
duj 1/2 - /off(u;) 



27r [u: - eo)2 + F2/4 
X T{u:)[f{u:-^iL)-f{uJ-^iR)], (28) 

/off(w) = E Tr /(t^ - Air)/r, 
r=L,R 

are consistent with those calculated by Hershfield [Eqs. 
(87) and (88) in Ref.^^], except for a factor 2 and an 
offset. The factor 2 is lost since we consider the spinless 
case. The offset would be related with the shift appeared 
already in the average Eq. pTjl . 



B. Exponential distribution 

In the limit of ig ~* oo, the number distribution follows 
from the inverse Fourier transform of the characteristic 
function Z within the saddle point approximation: 

Solid lines in Fig. [3] (a) are the equilibrium probability 
distributions of electron number in the weak tunneling 
regime derived from Eq. (j22p for the symmetric coupling 
(el/ = and Fl = Tr). When the dot-level is at the 



lead chemical potentials, eo = (line A), the distribution 
is well fitted by a Gaussian distribution around N Ki 
(dashed hne), 



(a) 



P{N) w exp 



-t. 



iN^{{ST))/to 



l21 



2{{Sr^))/to 



(29) 



where the first and second cumulants are given by 
Eqs. dM]) and (gS]). Around iV = ±1/2, it is suppressed 
strongly. Since our RLM possesses the particle-hole sym- 
metry, the distribution appears symmetric for negative 
and positive eo (not shown). 

The deviation from the Gaussian distribution grows 
rapidly as the dot-level leaves away from the chemical 
potential. The range where the Gaussian distribution is 
valid shrinks (lines B and C). Especially, for =Feo ^ T, 
the exponential- like distribution appears (line C); 

P(7V)«exp[<ors(±Af-l/2)]. 

The exponent is the decay rate of an excited state, 
namely for =Feo S> T, the decay rate of an empty/occupied 
state. The strong asymmetric shape of the exponential 
distribution is in clear contrast with the distribution of 
the open chaotic cavity^^, which should be symmetric 
around N = (N) in equilibrium. 

In the limit of the occupied or empty QD, a more 
careful treatment on the branch cut of the CGF is re- 
quired. The square root in Eq. ((22l) reads \/D w 

y(±T~MC7rsF for T^a/T -^ oo. Then the branch 
is chosen uniquely from the normalization condition 
W(i'(0,0) = 0. We obtain W(i) = ±ito^/2 and thus 
the delta distribution, 

P{N) =(5(iVTl/2). 

It supports our intuition that without the thermal and 
quantum fluctuations, an electron (a hole) is localized in 
the QD. 



III. SPIN EFFECT IN AN INTERACTING QD 
WITHIN THE MARKOV APPROXIMATION 

Up to now, we omitted the spin degrees of freedom. 
The trivial extension of the Hamiltonian ([5]) is to intro- 
duce the spin index a =t, i for fermions as, ark — > CLrka- 
and d ^ da- The extension results in the CGF twice 
as large as Eq. ITOl) because of the spin degeneracy. A 
nontrivial spin effect appears when the on-site Coulomb 
interaction Hint = Udld^d]di is accounted for. For sim- 
plicity, we will consider the limit U—^oo, where the dou- 
ble occupancy is forbidden and thus the QD states are 
limited to | f), | |) and |0). Furthermore, we will limit 
ourselves to the weak tunneling case T ^ F, where the 
experiments have been done^i^ and the master equation 
approach^'^ works. 

In this section, first we construct the FCS theory for 
the joint probability distribution of current and electron 
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FIG. 3: The logarithm of the probability distribution of elec- 
tron number P{N) for the noninteracting (a) and interacting 
(b) QDs in equilibrium for the symmetric case (Vl = Fij). 
The vertical axes are normalized by the coupling strength F 
and the horizontal axes are shifted by —1/2. Solid lines are 
for various dot-levels. Dashed lines are for Gaussian approxi- 
mations. 



number within the Markov approximation following the 
theory by Bagrets and Nazarov^'- fSec. IIIl!^ . Then we 
will discuss the spin effect in Sec. IIIIBI 



A. Extension of the FCS theory in the master 
equation approach 

Within the Markov approximation, a particular time 
evolution of the QD state is characterized by s tunneling 
events happen via a junction r^ into/out of {vi — ±1) the 
QD at successive times in a time interval (— ioo < ^i < 
t2 < ■ ■ ■ < ts < toe)- Such a sequence of the events, the 
sample, is then written as, 

Cs = iti,ri,vi;t2,r2,V2;- ■ ■ ;ts,rs,Vs). 

In the following, we assume too^to^l/T. The sample 
t^s determines the time evolution of the QD state Iq -^ 
li ■ ■ ■ Is-i — > /s- It also determines the time evolutions of 
the current through the junction r into the QD and the 
electron number inside the QD as. 



Ir{t,Cs) = 


= eJ^^Vj S{t-t j)dr,r,, 


(30) 


nitXs) '- 


S 


(31) 



Here Ni is the number of electrons inside the QD for the 
state I. The probabihty to find the QD in a state m 
evolves according to the master equation, 



dtPm{t) = -y^^LrnnPnjt). 



The off-diagonal components of the transition matrix L 
are given by the tunneling rate within Fermi's golden rule 
Eq. dSSl), 



Lmn — — 2^ ^r ^N^,N„±1- (32) 

r=L,R 

The diagonal components satisfy, 

Lnn = - 2J ^-'nn =l{n)- 



(33) 



The probability for the sample Cs to occur 
Qs(XU,ri,Vi\) = Q{ti,ri,vi]- ■ ■ ;ts,rs,Vs) is then 
determined as, 

go = e-270o)*»^ 



X e-'^('«)(*i+*°°), (34) 

where ts+i=tao- In the Markov approximation, the cur- 
rent I{t) and the number n{t), Eqs. H]) and 0, are 
simply numbers and commute at different times. Then 
the characteristic function of the joint probability distri- 

bution of qr = J_^/ i^ dt Ir{t)/e and r is, 

Z(A,0 = (O)j,, 0^explj2'^rqr + tCTy (35) 
Here the average is performed over all samples, 
(O)^ = OoQo+^ fdtiO{ti,n,vi)Q{h,n,vi) 



ri.vi 
toe. rt2 



+ y^ dt2 dti02{{U,ri,Vi})Q2{{U,ri,Vi}) 



{ri ,Vi} 



For the sample Cs, the stochastic variables Qr and r read, 
qrsi{U,ri,Vi}) = Y,f^J^^.^v-jSr,r, and Ts{{ti,ri,Vi}) = 
E'°7„.„ Ni^ {tj+i -tj) + 0(l/r) from Eqs. §0\) and (01]) 
and jmin(max) IS the minimum (maximum) j satisfying 
the condition \tj\ < to/2. The stochastic variable O then 
reads, 

Os{{ti,ri,Vi}) 
= Y]_ "^xp X! * '^'' ^J ^'■''■j + * ^ ^'j (*i+i ~ ^j) ' (36) 

i=jmin \ I- / 



for io > 1/r. By comparing Eq. ^ and Eq. ([361), we 
observe that the characteristic function ((35|) can be ex- 
pressed in a simple form as. 



(0>,, = Qo'*+5] / diiQ''«(ti,ri,^i) 

r-too /•t2 

dt2 / 

-ioo "'-to 



E 



dti Q2'^{{U,ri,Vi}) + 



(37) 



where Qs'^ is obtained from Qs with the following re- 
placement for jmin < j < jmax: 

p., ^j.„,^..,A.^^ 7(^,)-7(/,)+*^iV,,. (38) 

The first replacement was pointed out by Bagrets and 
NazaroAiii^: The counting field for the current can be 
absorbed in the off-diagonal components. Our result is 
the second expression: The counting field for the electron 
number can be absorbed in the diagonal components. 

To proceed the calculation, the formal similarity be- 
tween the master equation and the Schrodinger equation 
is utilizedii^: dt \p{t))^-L \p{t)), where {n\p{t)) = p„{t) 
and (n|L|TO) = Lnm- The operator L can be decomposed 
into three parts corresponding to Eqs. ((32|) and p3|) . 



L = j- 



E(r+ + f;). 



(39) 



Then, Eq. ([37|l is rewritten as the Dyson series, 



2i\,o =E<™i^^^p(~ / "'^^^w 



\nn) 



« ^(TO|exp(-ioLA.e) bo), (40) 

m 

where L{t) — L\^^ for \t\ < to/2 and L otherwise. To 
obtain the second line, we utilized the assumption to 3> 
I/r and the condition X)m("^l L = 0. We also assumed 
that the stationary state | po) , namely the zero eigen state 
L\pq) = 0, can be reached well before t = —to/2. The 
operator L\(^ is obtained from L with the replacement 



^A,^ =1 + 1 



e5]A^H0(/|-E r^^'^'^ + f^e 



T— „— i Xj 



(41) 
In the limit of to ^'^ oo, the minimum eigenvalue A,„in of 
the operator La,^ is dominant -Z«exp(— to Amm)- Thus 
the CGF reads, ' 



W(A,e)--toA„ 



(42) 



It is noticed that the normalization condition W{0, 0) = 
is satisfied since Amjn is expected to be the zero eigen 
value for Ar = C = 0. 



Let us apply the above mentioned scheme to the RLM. 
The QD states are hmited to the empty |0) and singly 
occupied |1) states. Then, \p{t)) is given by *(pi(i),po(0) 
and the operator (|4ip reads, 



iA.C = 



r±(A) 






p±g±<Ai 



-r+(A) 
r^e±*^«, 



(43) 



in the matrix form. Here we shifted the origin [Ni in 
Eq. (|¥T|) is replaced with Ni ~ 1/2]. As already demon- 
strated, unlike the counting field A in the off-diagonal 
components, the counting field ^ appears in the diagonal 
components: =F i ^/2 for the occupied/empty state. Two 
eigenvalues of the matrix (P5|) are 

A±-rs(l±\/i^)/2. 
From A_ and Eq. (|42|) . we reproduce Eq. ([22]) . 



B. Spin effect for U ^ oo 

Let us discuss the spinful case in the presence of strong 
Coulomb interaction based on the scheme in the preced- 
ing section. The probabilities to find the available three 
states, It), li) and |0) are given by *(pt(i),p; (t),po(t)) ■ 
The corresponding transition matrix with the counting 
fields reads, 



--i^/2 -r+(A) 

^A.c=| r--i^/2 -r+(A) 

-r-(A) -r-(A) 2T++i^/2 



(44) 



We found a unique eigen value Amin satisfying the 
normalization condition. The resulting CGF becomes 
Eq. dH]) with r+ replaced with 2r+. 

Figure [3] (b) shows the probability distribution of elec- 
tron number in equilibrium. We observe the exponential 
distribution both for eg^T and —eg^T (lines E and 
C), but with different exponents by a factor 2: In the 
former case, the excited state is 1 1) or 1 1). Therefore the 
decay process is the outgoing process of the spin inside 
the QD and the decay rate is rs«r. In the latter case, 
the excited state is the empty state |0). Thus the decay 
process is the incoming process of an up- or down-spin. 
Then the decay rate is the sum of up- and down-spin 
tunneling rates rx;~2r. 



IV. JOINT PROBABILITY DISTRIBUTION OF 
CURRENT AND ELECTRON NUMBER 

In this section, we discuss the joint probability distri- 
bution P{I, N) in the limits of large source-drain voltage 
and weak tunneling (/i^ = —^J'R and eo = for F ^ eV 
at T = 0). This case, the dot-level is inbetween the 



two chemical potentials, and Eq. (|22p is also valid even 
at zero temperature. These conditions were realized in 
the experiment^, but only the current distribution was 
measured. Figures S] (a) and (b) are contour plots of 
lnP(/, N). The inverse Fourier transform of the charac- 
teristic function Z is performed within the saddle point 
approximation, 

lnP(/,7V) « >V(i)(A*,r) -italX* -HqNC, 

For the symmetric case (a), we observe a clear feature: 
For /wO, the distribution of N is almost uniform, while 
for / ^ (/) a peak appears around N — 0. The for- 
mer corresponds to the rare tunneling case. In this limit, 
only one tunneling event occurs during to and one kink 
appears in the experimental data [upper panel of Fig. [Tj . 
Since the position of the kink is arbitrary, r takes any 
value between —1^/2 to io/2 with the same probabil- 
ity. The latter corresponds to the frequent tunneling case 
[lower panel] , where many kinks are randomly distributed 
and thus r w to{N) = 0. We would like to stress that it 
is not obvious that the crossover between the two be- 
haviors occurs at / = (/) and N — (N) as demonstrated 
in Fig. [4] (a). For the asymmetric case, the above de- 
scribed features are modified [Fig. [4] (b)] . Equidistance 
contours around / w implies that the electron number 
is exponentially distributed. At the same time, a long 
tail appears in the regime />(/), which implies that the 
current is Poisson distributed. 



EFFECT OF LARGE QUANTUM 
FLUCTUATIONS 



For now, we have mainly analyzed the number dis- 
tributions and the joint probability distributions within 
the Markov approximation in Sees. Ill 51 IIIII and [TVl In 
this section, we turn our attention to the non-Markovian 
effect enhanced by increasing strength of tunneling. Sec- 
tion |V3 sketches the derivation of the CGF for the SET. 
In order to go beyond the Markov approximation, the 
Keldysh generating functional approach is adopted again. 
Then in Sec. lVBl we will compare the probability distri- 
butions of electron number for two cases, the RLM and 
the SET. 



A. SET in the intermediate conductance regime 

Here we consider the SET, the metallic QD. The energy 
levels inside the metallic island are continuous. However, 
the single-electron charging energy Ec can be larger than 
the temperature and thus the low-energy physics is dom- 
inated by two island states, the charge neutral state |0) 
and the state with an excess charge e |1). The energy 
difference between the two states Aq = Ec{l — 2QQ/e) 
is controlled by the gate-induced charge Qq (Aq corre- 
sponds to eo for the RLM). The Hamiltonian is written 




FIG. 4: Contour plots of the joint probability distribution 
lnP(7,A'') are given for the symmetric case in (a) and the 
asymmetric case in (b) {Tl —bFit). The function lnP(/, A'^) 
has a maximum value In P{I, N) = at I — {1} and A'' — (N) = 
(r_L— r_R)/2r. The counter interval is for/20 and the shaded 
regions are for very small values. Vertical axes are shifted by 
-1/2. 



with a pseudo-spin-1/2 operator acting on the two states 
a+ = |l)(0|and<7, = |l)(l|-|0)(0|as2a, 



r—L^R,I km 

+ Yl J2^^^^lkm"'rk'rn^++ii-C-). (45) 

r^L^R kk' m 

The operator a]!/,„ creates an electron with the wave vec- 
tor k in the left or right electrode or the island (r=L,R,I). 
The transverse coupling describes the tunneling, which is 
assumed to preserve the spin and the transverse channel 
(The channels are denoted by m). The junction conduc- 
tance is written with the tunneling matrix clement T^, 
the number of channels A'ch and the electron DOS pr as 
1/Rr = 2^6^ NchlTrl"^ piPr- A measurc for the tunneling 
strength is the dimensionless conductance ao = a^ +ctQ, 
where Oq = l/2TTe^R,, 
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are introduced following Refi^: 



h{t). 



(46) 



They are fixed in the same way as for the RLM case, 
Eqs. (O and ([5]). The voltage drop in the junction 
r, /ir, is determined by the capacitances of the left 
and right junctions Cl and Cr as Pl{r) = '*l(_r)C^ 

[KL(fl)=±Cfl/i/(CL-HCfl)]. 

The following calculations are completely in parallel to 
those in Ref32., therefore we will just briefly sketch them. 
The starting point is the effective action of the pseudo- 
spin for TVch —^ oo derived from the Hamiltonian (j45p 
after tracing out the electron degrees of freedom. By 
introducing the Majorana representation's,, az = 2ffc— 1 



and (Tj 



St. 



1, where c and 4> are Dirac and Majorana 



fermions, the parts describing the charging energy ^ch 
and the tunneling St are written a o^°i^^'^'^ , 



dt c{ty [idt - Ao + h{t)] c{t) + -cj){t)dtm 



Sc\i 



St = - dtdt' c*{t)(t){t)a{t,t')(t>{t')c{t'), 
Jc 

where c and (j) are the complex and the real Grass- 
mann fields, respectively. It is the Bose Kondo model 
with anisotropic coupling and the magnetic field. The 
bosonic propagator, the particle-hole GF a(f, t') — 
J2r=L,R'^rit,t')e''^'^''^*'>-'^-^*">^ dcscribes the tunneling of 
an electron from the lead r into the island. The bare part 
is Ohmic except for the offset and written in the Keldysh 
space as, 



arito)' 



-ITT a. 



° {u-prY+El 




The high-energy cutoff Ec is introduced to regularize 
the UV divergence. The path-integral is performed by 
introducing a Grassmann source field J, 

Z = D[c*,c,<j>] exp{iSch + iSt) 



exp 



E 



(-1)'^ 



Tr 



S 6 



X expl - / dld2J(l)g^(l,2)J(2) 



»Tr[lns7'] 



,7=0 



The trace is understood as time integration over C . The 
Dirac and Majorana GFs, gc and g^, satisfy, 

[idt'A(> + h{t)]g,{t,t') = 5(i,t'), 
idtg4t,t') = 5{t,t'), 

and the anti-periodic boundary condition such as 
Eq. pi)) . Further calculations are based on the diagram- 
matic expansioii^^. The diagrams in Fig. [5] are resummed 
in our approximation^Si^i^^. 



W 



RTA 



The source fields for current and number, 



Vr 



and ft.. 



Trln(iG-^) 



G-^ = g-' - s, (47) 

~ig^{t',t)a(t,t'). (48) 



We can see Eq. p7|) is formally the same as Eq. (TTSl) . 
Then it is almost clear that the result is formally the same 
as Eq. (fT9|l , but the Keldysh and retarded components of 
the self-energy corresponding to Eqs. P^ and p?]) are 
replaced with those of the SET: 



Ef(u;) . 


-- 2aS 






vrT 




- zTrao (w — ^r)coth 
V 


2T J ' 




Sf (a;) = 


= -2i7r 


al (w - ^^), 







where ?/; is the digamma function. Consequently, the 
transmission probability Eq. (|^D|) reads, 



T(c.) 



(271")^ nr=L,fl,"o('^ - Mr) COth 



to — ^r 

2T 



Ao~Er^L,R^H^) 



It reproduces the transmission probability derived in 
Rcf.''^, which can describe the higher order inelastic co- 
tunneling'^— as well as the sequential tunneling. 



o 




©=o. 



FIG. 5: Diagrammatic expansion of W. Solid, dashed and 
wavy lines are for Dirac, Majorana and particle-hole Green 
functions. 



B. Probability distributions of electron number for 
SET and RLM 

In the hmit of ao^O, the CGF W^^^ reproduces that 
of the extended master equation approach. The formal 
expression is the same as Eq. (22) with the tunneling rate 
of the RLM, Vf, replaced with that of the SET, 



rl I Ir 



±- 



Ao-/ir 



^i?,. (e 



±(An- 



-t^r)/T_iy 



Therefore within the Markov approximation, though this 
time the tunneling rate depends on the source-drain volt- 
age, the current and number obey qualitatively the same 
statistical distribution as that of the RLM. 

However, it is not the case when the tunnel coupling 
increases and the quantum fluctuations are enhanced. 
We compare the nonequilibrium number distributions for 
the RLM [Fig. [6] (a)] and for the SET [Fig. [6] (b)] at 
eo = Ao = in the symmetric case (Vl=Tr, ct'L = ct% and 
fJ.L — —fJ'R)- At very small tunnel coupling, the distribu- 
tions for both cases approach the master equation results 
[dashed lines overlapping with curves for r/el^ = 0.05 in 
panel (a) and for ag = 10~^ in panel (b)]. Here, vertical 
axes are normalized by Fs at OK, namely F for the RLM 
and {Rl + RR)V/2eRLRR for the SET. For the large 



tunnel coupling {T/eV and ao), the fluctuations of the 
electron number get stronger and thus the distributions 
for both cases shrink. However, a qualitative difference 
appears: For the SET, the width of the distribution just 
shrinks, but for the RLM, a two sided exponential dis- 
tribution like profile shows up. We can expect that it is 
related with the difference in the dominant higher order 
tunneling processes. For the RLM, the large tunnel cou- 
pling enhances the multiple tunneling process of a single 
electron, which preserves the phase coherence. For the 
SET with large Nch, the coherent process, the elastic co- 
tunneling process, is suppressed as compared with the 
inelastic cotunnehng process^. In the latter process, dif- 
ferent electrons tunnel back and forth between the leads 
and the metallic island. Therefore, the phase coherence 
is lost. 
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FIG. 6: The nonequilibrium number distributions are given 
for the resonant-level model in (a) and the single-electron 
transistor in (b) at eo = Aq = and T = 0. Solid lines are 
for various coupling strengths F in (a) and for various dimen- 
sionless tunnel conductances ao(= Z^r=L r l/27re^/?r) in (b). 
The dashed lines overlapping with curves for V jeV = 0.05 
(a) and for ao = lO""* (b) are obtained by the extension of 
the master equation approach. Horizontal axes are shifted by 
-1/2. 



The phase coherence effect is pronounced when the 
QD conductance approaches to the quantum conduc- 
tance e^/(27r?i). In this regime, the high band- width 
of the QPC electrometer is required - even for a small 
source-drain voltage, 0.1 mV, the band- width is required 
to be over than 10 GHz. For accurate measurements, the 
QPC current should be much larger than 1 nA, which 
would make the back action neglected here a crucial in- 
gredient'^^. The FCS theory of electron number including 
the QPC electrometer will be discussed elsewhere. 
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VI. SUMMARY 

In conclusion, we have extended the FCS theory for the 
statistical distribution of electron number inside the QD. 
We found the non-Gaussian distribution, the exponential 
distribution, when the dot-level is far away from the lead 
chemical potentials. The exponent is proportional to the 
spin degeneracy. The large tunnel coupling would nar- 
row the number distribution generally. Especially for the 
symmetric noninteracting QD, the two sided exponen- 
tial distribution appears in the probability distribution 
of electron number out of equilibrium. One measurable 
prediction for the presently available experimental setups 
is the joint probability distribution of current and elec- 
tron number, which reveals the nontrivial correlations 
between the two quantities. 

In the present paper, we have extended the FCS the- 
ory of the master equation approach^'^ and derived the 
general expression for the CGF of the joint probability 
distribution. We have shown that the CGF is obtained 
from the minimum eigenvalue of the transition matrix 



with the counting field for electron number in the diago- 
nal components (Sec. IIIlX| . The non-Markovian effects 
have been treated in the frame of the Keldysh generating 
functional for the two cases, the non-interacting (RLM) 
and the interacting (SET) QDs. For the RLM, it is pos- 
sible to derive the exact expression (Sec. Ill A|) . For the 
SET, we performed resummation of diagram a^°i^^'^'^ to 
obtain an approximate CGF in the intermediate conduc- 
tance regime (Sec. IV A[) . 
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